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We derive a generic phenomenological model of a Majorana Josephson junction that accounts for 
avoided crossing of Andreev states, and investigate its dynamics at constant bias voltage to reveal 
an unexpected pattern of any-7r Josephson effect in the limit of slow decoherence. 
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Recently, the proposals of solid-state realizations of 
Majorana fermions came into focus of attention. While 
the first proposal [l[ concerned non-Abelian excitations 
in 5/2 FQHE in semiconductor heterostructures, most 
proposals 0, Q exploited exotic superconductors where 
Majorana fermions correspond to zero-energy states of an 
effective BdG Hamiltonian. The Majorana states are in- 
strumental for realization of topological quantum compu- 
tation Q • More recent contributions [B[ utilize the prox- 
imity effect from a conventional superconductor, either in 
nanowires in a strong magnetic field and with strong spin- 
orbit interaction [6|-|8| , or in topological insulators [{J [lOj . 
This brings the Majoranas close to experimental realiza- 
tion, and underlines the importance of reliable experi- 
mental signatures of their presence. Among the signa- 
tures are half-integer conductance quantization llj and 
47r Josephson effect in superconductor-superconductor 
(SS) junctions0,H[l3|. 

No 4-7T periodicity is to be seen in the stationary ground 
state of the junction It can only be observed [lj, [l5[ in 
dynamics induced, for instance, by a d.c. voltage bias. 
Unambiguous signature of this anomalous periodicity is 
the current noise peak at half of the Josephson frequency 
Uj = 2eV/H We have suggested that the avoided 
crossing of Andreev states is intrinsic for finite systems 
and restores the 2ir periodicity of the junction ground 
state |16j. This has been recently confirmed by detailed 
calculations of the Andreev spectrum of the nanowire- 
based SS junctions [171 ] . 

In this Letter, we put forward a generic phenomenolog- 
ical model of a Majorana Josephson junction and demon- 
strate that the dynamics in the junction are substantially 
richer than thought. In particular, the sharp peaks in 
noise spectrum of a voltage-biased junction are not gen- 
erally confined to any definite fraction of ujj : one can talk 
of any — it Josephson effect in this context. Experimen- 
tal observation of these singularities would give a robust 
proof of the existence of Majoranas and open up the pos- 
sibilities for quantum manipulation of these states. Our 
treatment of dynamics encompasses the Landau- Zener 
tunneling at the avoided crossings, decoherence, relax- 
ation, and quasiparticle poisoning. 

We exemplify with a nanowire setup (Fig. [IJ although 
the same phenomenology extends to topological insula- 




FIG. 1: a. A Majorana Josephson junction is formed by 
mounting a nanowire (black) on two superconducting leads 
(grey) resulting in four Majoranas 71-4. b. The energies 
of the junction states versus phase near an avoided crossing 
point, c. Corresponding Andreev levels, d. The energies at 
bigger phase scale. 



tors. A nanowire mounted on a single superconducting 
lead develops a topologically non-trivial state in a pa- 
rameter range of magnetic fields and gate voltages [7]. 
Two Majorana states emerge at the wire ends. Majo- 
rana Josephson junction is formed by mounting the wire 
on two leads biased with superconducting phase differ- 
ence </>. Two extra Majorana states 72,3 emerge at the 
junction, in addition to the end states 71,4. The overlap 
between 72 and 73 is strong but does depend on phase 
and vanishes at a certain <f>Q. If one disregards the end 
states @,[l2, 13], the resulting energies are An periodic in 
(j> and the resulting states are of indefinite parity. We ex- 
emplify this dependence with E(<f>) = ±e sin( ^~^ ). To 



fix the parity, it is paramount to bring the end states to 
the picture. We developed [l6| a scattering matrix theory 
where the 2tt periodicity is proven from the topological 
properties of the scattering matrix. In a nutshell, the 
crossing of Andreev levels is avoided. We need a practi- 
cal Hamiltonian to describe the details of the situation 
in the vicinity of 4>o- That can be rigorously derived 
from the scattering approach, yet we opt here for a sim- 



pie heuristic deviation in terms of overlaps of Majorana 
states. 

These overlaps are exponentially small for long wires, 
oc exp(— L/2£), L being the wire length, the localiza- 
tion length £ being of the order of the spin-orbit length 
L so . For InAs wires [18|, L so = 0.2/im, and L would not 
exceed 2\im since inevitable disorder forbids topological 
state for longer wires. This sets the biggest exponential 
suppression to ~ 1CP 2 . Owing to the exponential sup- 
pression, the direct overlap ii4 between the end states is 
much smaller than the overlaps between the end and the 
junction states, and can be disregarded. This brings us 
to the following Hamiltonian: 



H = e((j>- (j) )j2j3+ 

(tl272 + il373)7l + (*42 72 + ^4373)74- 



(1) 



that is valid in the vicinity of the crossing point and pro- 
vides a generic phenomenological model of a Majorana 
Josephson junction. Here, the overlaps t are real, and 
7i_4 are self-conjugated anticommuting Majorana oper- 
ators [ijj. 

It is instructive to give the eigenenergies of the full 
many-body states of the Hamiltonian, rather than the 
associated Andreev levels. The Hamiltonian conserves 
the parity of the particle number and therefore gives rise 
to eigenstates with either odd or even number of parti- 
cles. There are two eigenvalues of opposite sign for each 
parity, 



±\/G» e + -( 



where 



Go,e — 5 v (tl2 ± t 



43 ) 



(*13 T U2Y 



(2) 



(3) 



and ± sign is chosen such that G e > G . Their phase 
dependence (Fig. [TJd) gives a familiar glimpse of avoided 
level-crossing hyperbolas, G o . being the minimum en- 
ergy splittings of odd/even states, respectively. The two 
positive energies of the associated Andreev levels are 
given by £1,2 = \E e \ ± \E a \ (Figflb). This characteristic 
form is conformed by the numerical calculations based 
on microscopic models proving the validity of the 

Hamiltonian ((TJ) . The absence of the direct overlap tu 
leads to a special property: the phase-dependent term 
describing the overlap of 72, 73 anticommutes with the 
rest of the terms. This guarantees the energies to be 
even in phase and to merge far from <p = <pQ, these prop- 
erties would be absent for a most general four-Majorana 
Hamiltonian. 

Let us notice that the junction in either odd or even 
state is nothing but a qubit that is similar to other super- 
conducting qubits that commonly exploit avoided level 
crossing One can employ quantum manipulation 

of the resulting Majorana states by changing the super- 
conducting phase in time. For instance, following [20j, 
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FIG. 2: Processes affecting the junction dynamics. In between 
the parity-dependent Landau-Zener scatterings (described by 
2x2 matrices S , e ) the junction is subject to dephasing (with 
a rate T d ), relaxation (17) and quasiparticle poisoning (T q ). 



one can prepare the qubit in the ground state reason- 
ably far from the crossing point (f> — 4>i , and give a pulse 
that brings the junction to <p = <pQ. This will cause Rabi 
oscillations with frequency 2G/h that can be detected 
by measuring the probabilities to find the qubit in the 
ground or excited state after the pulse as functions of 
pulse duration. 

Here we restrict ourselves to the case of immediate 
experimental relevance where the junction is biased by 
a d.c. voltage V so that the phase <j> is swept linearly 
with time, 4> — 2eV/h. In a usual Josephson junction 
where the energy levels are 2n periodic, such bias re- 
sults in coherent oscillations of the supercurrent /(</>) = 
2e/HdE(4>)/d4> with Josephson frequency ujj = 2eV/h 
[2J|. The idea behind the 4n Josephson effect (25| is an 
apparent 4-7T periodicity of energy levels in the limit of 
vanishing G, this suggests the oscillations at a half of 
Josephson frequency, I(t) — ±I m cos(u)jt/2), I m = ee/H. 
Albeit these oscillations cannot be coherent owing to ran- 
dom switching between the two branches ± of the energy 
spectrum. The signature of An periodicity is rather a 
sharp peak in the spectral density of the current noise 
141] . with the width of order of switching rate, and inte- 
grated intensity being given by i^/2, the average square 
of the current. For this simplified picture to hold, one 
should require sufficiently small voltages, V -C e. Failure 
to satisfy this condition results in proliferation to higher 
energy levels and finally to continuous spectrum, this in- 
creasing the peak width to the values of the order of A 
and thus rendering noise peaks undetectable. [3, [ItJ 

In this Letter, we address the noise in Majorana 
Josephson junctions at smaller voltages. Evidently, the 
avoided level crossing results in usual Josephson effect in 
the limit V — > 0. The complex and interesting crossover 
between 27T and 47r regimes involves Landau-Zener (LZ) 
tunneling upon crossing a point <f> = <f>Q+2im in the vicin- 
ity of the point. The parity obviously does not change, 



3 



and for each parity we have a classic setup of LZ tunnel- 
ing [23| between two levels. The values of the qubit wave 
function before and after LZ scattering are related by a 
2x2 unitary matrix: 



where the probability of LZ tunneling is given by 



exp 



4tt Gl 
~ eV~ 



(4) 



(5) 



This suggests an importance of a voltage scale eVo =— 
AnG 2 1 1 <C e at which the probabilities are of the order 
of 1 and the crossover between 27r and 47r regimes is ex- 
pected. We stress that the probabilities are generally 
different for odd and even states that permits the iden- 
tification of these states that are hardly distinguishable 
otherwise. 

The quantum dynamics are affected by the processes of 
relaxation, dephasing and quasiparticle poisoning (Figj2]) 
that occur throughout the time-line with no peculiarities 
near the crossing points. We assume low temperature 
ksT <C e, so that the relaxation is always from higher to 
lower energy state with the rate F r ((f)). The decoherence 
suppresses the non-diagonal elements of the density ma- 
trix (with the rate Td((j>)) not affecting the diagonal ones. 
We assume the fluctuation of the phase <fi to be the main 
source of the decoherence, in this case Td(4>) oc I 2 (4>)- 
The quasiparticle transfer processes account for a parity 
change. They may be due to stray quasiparticles in the 
bulk superconductor that come to the junction with the 
energies of the order of the superconducting energy gap 
A > e and lose this energy either adding or annihilating a 
quasiparticle in Andreev levels under consideration. Due 
to significant initial quasiparticle energy, the probabili- 
ties to find the junction in either upper or lower state 
after a quasiparticle transfer, are the same. The quasi- 
particle rate T q does not depend on the phase 4>. 

This results in the straightforward but lengthy equa- 
tion for the density matrices p _d to be found in [26[. We 
solve this equation with continuity conditions p .d(t a c + 
0) = SoApoAtac - 0)5*^^, t ac corresponding to time mo- 
ments of the crossings, and compute the correlator of 
current operators 



dh P dhe^^Hilih)!^))). (6) 



S(u) 



that gives the spectral density of the current noise. We 
concentrate on two limiting cases of fast (Fig. [3]) and 
slow (Fig. [3J decoherence. In both cases, we assume slow 
relaxation and poisoning, ujj ^> r r ,r 9 . 

"Fast" implies the quantum coherence is lost during a 
period of the Josephson oscillations, Yd ^> Wj, and the 
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FIG. 3: The spectral intensity S(u>) of the current noise 
for a set of V corresponding to LZ probabilities shown, in 
the limit of fast decoherence. a. Indistinguishable parities 
P e = P = P LZ ; b. P c P (G c = 4G ). Distinct peaks at 
multiples of ojj at P 1 transmute into a single peak at w,/2 



at P fts 1. (r r 



0.02ojj for all plots). 



equation for density matrix reduces to a master equa- 
tion. Fig. [3^ shows the spectral density for equal LZ 
probability for even and odd states, G e ~ G D . The volt- 
age growth from the lowermost to upper curve resulting 
in increased Plz- At low voltage (Plz <§C 0, the noise 
peaks at utj as well as at its multiples, the latter mani- 
festing non-sinusoidal I(<fi). This proves a usual period- 
icity. At higher voltage where Plz ~ 1 we see a single 
peak at uj/2 manifesting 4tt periodicity. In both limit- 
ing cases, the peak widths ~ T rt q P . Important feature is 
the absence of any distinguishable peaks at intermediate 
Plz- The reason is the LZ tunneling causing incoher- 
ent switching at almost any crossing point. The peaks 
acquire width ~ uij 3> Y r ^ qp and correspondingly reduce 
their height to the background level. 

In Fig. (3t> the LZ probabilities are very different at the 
crossover corresponding to G e /G = 4. Now one can dis- 
tinguish the peaks at both u>j/2 and uij in the crossover 
region, though they are reduced in height in comparison 
with the limiting cases. The explanation is the parity sep- 
aration in time domain. Since T q -C 0Jj , the parity persist 
over many periods between the random switches. While 
the junction is in even parity state, Plz ~ 1, and during 
this time interval the noise at uij/2 is generated. While 
the junction is in odd parity state, almost no LZ tunnel- 
ing takes place, and the noise is generated at Josephson 
frequency. The experimental observation of two peaks 
would thus prove the parity effect. One can also think 
of a more challenging observation where the noise can be 
resolved fast, that is, at a time-scale < T^ 1 . Such noise 
measurement will monitor the parity of the junction in 
real time. 

The results in the opposite limit of slow decoherence 
Td *C uij are decisively more complex and intriguing (Fig. 
3]). In this limit, the dynamics are truly quantum over 



many periods. An analytical analysis gives the positions 
of the noise peaks as well as the integrated noise intensi- 
ties around each peak Most striking feature is an os- 
cillatory dependence of the peak intensities and positions 
on voltage. This is a manifestation of quantum inter- 
ference between the subsequent LZ tunneling events not 
suppressed by decoherence. Similar interference patterns 
have been predicted and observed for superconducting 
qubits in[22, We have found that a voltage-biased 
Majorana Josephson junction presents the simplest and 
most striking framework for this interference effect. 

The quantum phase 9 accumulated between the sub- 
sequent crossing points is estimated as 



S/So 



6= J dt 

period 



AE(cf>(t)) 



(7) 



The phase is big on the scale eV/e, its increment by 
2ir gives an estimate of the oscillation period in voltage 
AV = (7r/8)V(eV/e) < V. 

Importantly, the frequency positions of the additional 
noise peaks (Fig. [4^), those the main Josephson peaks 
at multiples of uij, are not at any integer fractions of 
ujj. In the context, we can dub this any-7r Josephson 
effect. It stems from a quasi-energy splitting in a pe- 
riodically driven qubit. At V 3> Vb, additional peaks 
converge at (2n + l)ujj/2 oscillation around this fre- 
quency. The spread of these oscillations Auj does not 
vanish with increasing V: rather, it increases following 
Auj ~ (2e/h)y/VVo- This proves that any-7r Josephson 
effect can be observed at voltages V 3> Vb far beyond 
the crossover region. The width of the peaks is deter- 
mined by IV From this, we estimate the minimum de- 
coherence rate permitting the resolution of the peaks: 
Td — (e/h)-\/VVo- For the sake of simple drawing, we 
assumed indistinguishable parities such that P Q = P c . 
If P a P c , the additional peaks split once again corre- 
sponding to the two parities. [26j 

At V <C Vb, the noise intensity is mainly concentrated 
at a main peak at u>j. In the opposite limit, the intensity 
concentrates at the peaks converging to uij/2 retaining 
oscillating features even at high voltage. (FigJ3b,c) 

To summarize, we have derived a generic phenomeno- 
logical Hamiltonian to describe a Majorana Josephson 
junction with avoided Andreev level crossing, and inves- 
tigated its quantum dynamics at constant voltage bias 
with emphasis on noise signatures of the anomalous pe- 
riodicity. While in the fast decoherence regime the sig- 
natures follow an expected pattern, the interference of 
the subsequent LZ tunneling events results in a com- 
plex any-7r Josephson effect pattern in slow decoherence 
regime. The experimental observation of the effects pre- 
dicted will provide unambiguous signature of Majorana 
states in Josephson junction and open up the perspec- 
tives of quantum manipulation and parity measurements 
in such junctions. 
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FIG. 4: The slow decoherence limit. We chose epsilon/eVo = 
30 for all plots, (a). Frequency positions of the noise peaks 
versus V (ujo = (2e/ft)Vb). The main ones are at nuj 
while the positions of additional peaks oscillate converging 
at (n + l/2)uij. Only n = 0, 1 are shown, (b) Integrated noise 
intensity (in units of So = im/2) of the first two additional 
peaks, (c) The same for the n = 1 main peak. 
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APPENDIX 

Periodic continuation of the phenomenological Hamiltonian 

Far from the avoided crossing point (\(f> — <po\ 3> G/e <§C 2ir), the energies merge together, so that the lowest Andreev 
level is close to zero, E% ~ G 2 /e « G C e. At further increase of \<f> — 4>q\, the energies deviate from their linear 
asymptotes near the crossing points to become 2w periodic(Fig. [T}i). We approximate their dependence with 



E e = ±e 



(8) 



It may seem that the Hamiltonian may be extended to the full range of the phases simply by replacing tf> — cf>o with 
2sin((0 — <po)/2)). However, this would result in a 47r-periodic Hamiltonian, this being at odds with the natural 2tt 
periodicity. 

To resolve this apparent discrepancy, we note that the choice of Majorana operators in Eq. [1] that describes the 
lowest energy states, is not unique: one can substitute instead of 71-4 any linear combinations of Majorana operators 
of a bigger set that obey the commutation relations. The members of this bigger set would correspond to higher 
energy states not considered. The choice of four Majorana operators 71^4 made does not depend on phase in the 
vicinity of the avoided crossing point and does depend at a bigger scale. The choice in fact is a 47r periodic one. To 
give an example of a 27r-periodic Hamiltonian, let us substitute 73 in the form of two Majorana operators that do not 
depend on the superconducting phase, 

73 = cos((0 - M/2)% + sin((0 - 0o)/2)73 (9) 
With this, we can rewrite a seemingly 4tt periodic term in the form 

esin((0 - o )/2) 7 273 -> \ ((1 - cos(0 - + sin(0 - faYft) (10) 

which makes the 2tt periodicity explicit. In further considerations, we will work in a basis of the energy eigenstates 
that is explicitly 2tt periodic. This makes irrelevant the details of Majorana representation outlined here. 



Equation for density matrix 

Let us give here the evolution equations on density matrix that are straightforward but too bulky to fit into the 
main text. 

Beyond the vicinities of the crossing points (f> = (f>o+2im, n being an integer, we can disregard the terms proportional 
to G ,c- We work in local eigenenergy basis. To denote the elements of the density matrix, we use the subscripts e 
and o for even and odd parity sectors, respectively, and u and I for upper or lower states in each parity sector. We 
collect all the incoherent processes: relaxation, dephasing and quasiparticle poisoning. We introduce energy splitting 
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E((f>(t)) = 2? sin 



With this, the equations read: 

dpeu,eu 



dt 

dp el, el 

dt 

dpeu,el 

dt 

dpel,eu 

dt 

dpou,ou 

dt 

dp i,oi 
dt 

dpou,ol 

dt 
dt 



r rPeu.eu qPeu,eu \Pou,ou + Pol,ol) 

r V Peu,eu ^qPeu,eu 2^*2 i^Pou.ou Pol : ol) 
(-iE(4>(t)) -Td-Vg- |r r ) p eM , e i, 
(iE((j)(t)) - T d - T q - pei.eu, 
^rPo-UjOtt rgP ou ou -|- (Peu,eu H~~ Pel,el) 

PoUjOvX^r Poi^oii^g 2^'? (.Pgu,cu ~l~ PeZ,ez) 

(-iE{4>{t)) - T d -T q - |r r ) /0 ou , o ; 
(iE(4>(t)) - T d -T g - ir r ) p L 0tl . 



(11) 
(12) 

(13) 
(14) 
(15) 
(16) 
(17) 
(18) 



To specify the dependence of the decoherence rate, we assume that the decoherence takes place mainly due to the 
non-ideal bias conditions that give rise to the fluctuations of <f>. Those cause the fluctuations of the energy splitting 
SE = (dE /d4>)d(f> and are converted to the fluctuations of quantum phase thereby. In this case, T d oc (dE/dcj)) 2 and 
can be thus chosen to be of the form Td((f>) — r^cos 2 ((0 — 0o)/2). 

In the vicinities of the crossing points we may neglect the incoherent terms. For each parity, the Hamiltonian ([TJcan 
be written as a 2 x 2 matrix 



H = 



G„, c e(0 - 0o ) 
-l{4> - 4>o) -G Q ,e 



Integrating this over time and transforming to the eigenenergy basis reproduces the boundary condition on the density 
matrix given in the main text, for density matrices Pb(p a ) before (after) passing the crossing point, 



Pa — SpbS ; 

where the scattering matrix of LZ tunneling is given by Eq. U and x — § + Arg 
here is the gamma function) 



r(i 



+ 



,(ln 



eV e \ — eV e 



(19) 

i) (r 



Fast decoherence limit and master equation 

Let us first consider the case of the fast decoherence when the quantum coherence quenches at the time scale smaller 
than the period of Josephson oscillations, T d 3> ujj. This permits us to take into account only the diagonal elements 
of density matrix, the probabilities, describe their evolution with a master equation. We denote the probabilities with 
Pa, /3, where a = e, o at even or odd sectors and j3 = u, l for upper or lower state in the sector. With all the processes 
in Fig. [2] taken into account, the master equation for time intervals between the crossings reads: 

dPeU /-pi . -pi \ | -j i Pou + Pol /r, n \ 

—fa = -(IV + r q )Peu +T q - , (20) 



dpel 

dt 

dpou 
dt 

dpoi 
dt 



p, p, | t i Fou ^ Fol 

T r peu - T q p e i + T q - , (21) 

-(r r + r 9 K„ + r 9 ^i±^, (22) 

TrPou - TqPol + T q PeV ^ . (23) 
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This has to be supplemented by LZ boundary conditions at the crossing points for the probabilities before and after 
the passing. For the odd sector, we have 

j£l = P h ol + Po (Pla - P h ol) ; Plu = P h ou + Po {Pll - P h ou) ■ (24) 

Similar equation holds for the probabilities in the even sector. 

Let us label the four possible states ou, 61, eu, el with a single index j Solution of these equation in the long time limit 
approaches p^(t)j, that is periodic in time with the Josephson period. To compute the correlator of the currents, we 
also need the propagator of the evolution equation Uij{t 2 ,t\). It is defined at t 2 > h as the solution of the equation 
at the time moment t 2l Pi(t 2 ), with initial condition Pk(ti) = 5k j. 

The current is a function of a state given by 

h{t) = ImSgn(</> - <p ) cos((0 - o )/2) (<5i, ou + 5^ eu - S it0 i - 5 i)6 j) (25) 
with J m = (2e/H)(e/2), (j) — <fio + ujjt. The correlator is expressed as 

((/(i 1 )/(t 2 )))=^/ i (t 1 )/ J (t 2 )([/y(i 1) i 2 )-^(i 1 ))pf (t 2 ), (26) 
ij 

at t\ > t 2 , and is obtained by permutation of the time arguments otherwise. 

We solve the equation, the propagator and find the correlator numerically. The results for the current noise spectral 
density are presented in Fig. [3l 



Details of slow decoherence limit 



Interesting analytical results can be obtained in the opposite limit of slow decoherence such that Td,r,q <C u)j if in 
addition we assume that max (Plz, T q /ujj) 3> T r /ujj (the latter condition even in the absence of T q is satisfied at 
V > Vb/ln(r r /wj and thus certainly holds in the crossover regime). Under these assumptions, the relaxation can 
not set a preferential state. All possible states of the junction are present with equal probability, and long time limit 
density matrices do not depend on time and approach p D = p e = jl. In this case, we can neglect Td, r ,q implement 
the pure quantum dynamics to compute the current-current correlator at time separation |fi — 1 2 | <C r^J; . While 
not enough to resolve fine features of the noise spectral density such as the line-shapes of the noise peaks, this suffices 
to evaluate the integrated noise intensities in the vicinity of each peak. 

The positions of the peaks are not bound to the multiples or integer fractions of the Josephson frequency. To 
understand this in general, let us note that the independent solutions of the Schrodinger equation that is in our 
case periodic with the Josephson period Tj = 2ir/u)j, are Bloch-like functions of time satisfying 

|* fc (t + Tj)) = exp(iA fc )|* fc (t)) (27) 

and having Fourier components at discrete frequencies u}j(n + Xk/2n). The correlators thus can have Fourier compo- 
nents at all discrete frequencies satisfying uij(n + (A& — Xi)/2n). 

To analyze the quantum dynamics in the case under consideration it is proficient to apply a unitary transform that 
cancels the evolution of the wave function during the "free motion" between the crossing points. The phase difference 
X accumulated in the course of the free motion (Eq[7] ) is then ascribed to the LZ scattering matrix in certain parity 
sector so it becomes 

A _ ( Vl~Pe ie / 2 -e ix VP 




Since this matrix describes the evolution of the wave function over the period, its eigenvalues give exp(-iX) and 

cos(A) = VT~- Pcos(0/2). (29) 

Let us define 7 = arccos(-\/l — P cos(6/2))/ir, < 7 < 1, and recall that we have two parity sectors and correspond- 
ingly two parameters 7o, e - With this, the noise spectral density is digested in the form that makes the peak positions 
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explicit: 



S(ui) = (8((jJ — riLUj) + S(uj + nujj)) + 

n>0 

J2 S n°l W w - (7o + n)uj) + 5{u + ( 7o + n) Uj )) + 

n>0 

S n'°h i 6 i" - (~7o + n) Uj ) + S(uj + (- 7o + n)uj)) + 

n>0,± 

S£' e l (S(u - ( 7 c + n)wj) + S(uj + ( 7c + n)^)) + 
S n'°l i 6 ^ ~ (-7e + n)uj) + S(lj + (- 7c + n) Uj )) . 

n>0,± 

The integrated intensities S n of the main peaks are contributed by both parity sectors, while the additional peaks 
Sn '°^ c are in general resolved in parity with respect to the positions and height. It is worth noting that the definition 
of noise density in use is " quantum" , so that the noise at positive and negative frequencies does not have to be the 
same indicating the difference between emission and absorption of quanta with energy fvjj. However, this is not the 
case in the present framework: the calculation explicitly gives the spectral density that is even in frequency. 

To compute the intensities, we need to evaluate the correlator in Eq. [5J We note that in the representation used 



/(*) = I m cm(ujt/2)a z ; 



(30) 



(crossing points corresponding to t — nTj), and the evolution matrix £/(i2,ii) that gives the wave function at ti from 
the initial condition at t\ is simply given by 



f7(*2,*i) = {S) 



N 



(31) 



N being the number of the crossing points at the interval (t2,ti). 
For the intensities of the main peaks, this gives 



n _ jA 
J n — 1 m 



16n 2 



7T 2 (4n 2 - iy ^ 



sin 2 (6»/2) +P 0jC cos 2 (0/2) 



(32) 



As to the additional peaks, their intensities are given by 



4 cos(nD,/ujj)flLUj 



J o,e 

sin 2 (0/2) +P , C cos 2 (0/2) 



(33) 



where Q is the (parity-dependent) frequency position of the additional peak, P is the LZ probability at the corre- 
sponding parity. 

Making use of the above relations, we plot in Fig. U]the positions and intensities of the three lowest peaks assuming 
P Q = P e , that hinders parity resolution. 



Extra figures 

To illustrate the slow decoherence limit in more detail, we present here two extra figures. Fig. [5] illustrates the 
noise intensities of the two peaks converging to f2 = (3/2)to>j and of the second main peak at 51 = 2ujj. In both cases, 
the relative intensity is substantial in the crossover region V ~ Vq and slowly falls off upon increasing V. Fig. [B] 
represents the generic case of parity separation. The LZ probabilities here are taken to be distinctly different in the 
crossover region corresponding to G e — 4G . We see that a so-to-say 4-7T periodic peak at (Jj/2 in fact consists of the 
four distinct peaks of different intensity slowly converging to 1 /4 of the full intensity in the limit V 3> Vq. 
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FIG. 5: Intensities of the high-frequency noise peaks. The parameters are the same as in Fig. 2] Left: two additional peaks 
with frequencies converging to (3/2)a?j at V S> Vo .Right: main peak at 2wj. 




FIG. 6: The generic case of parity separation. Here, Vo = AtvGo/S, luq — 2eVo/h, e/cjo = 30, So = a. Peak positions. 

Four distinct noise peaks converge to half-integer multiples of Wj upon increasing voltage, b. Spectral intensity of the main 
peak at ojj. It is contributed by both parities, c. The same for the first two even parity peaks, d. The same for the two odd 
parity peaks. 



